PARTICLE APPROXIMATION OF SOME LANDAU EQUATIONS 
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Abstract. We consider a class of nonlinear partial-differential equations, including the spatially 
homogeneous Fokker-Planck-Landau equation for Maxwell (or pseudo-Maxwell) molecules. Con- 
tinuing the work of [6] [7] [4], we propose a probabilistic interpretation of such a P.D.E. in terms 
of a nonlinear stochastic differential equation driven by a standard Brownian motion. We derive 
a numerical scheme, based on a system of n particles driven by n Brownian motions, and study 
its rate of convergence. We finally deal with the possible extension of our numerical scheme to 
the case of the Landau equation for soft potentials, and give some numerical results. 
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1. Introduction and main results 

1.1. The equation. Let Sd be the set of symmetric c? x o? matrices with real entries, and S'^ its 
subset of nonnegative matrices. For a : M.'^ h- s- S'^ , we consider the partial differential equation 

(1) dtft{x) = ^ y] "^M / ""y^^ ~ y)\ ft{y)dj ft{x) ~ ft{x)djft{y)\dy 

where 9* = 9^ = and where the unknown {ft)t>Q is a family of probability density func- 
tions (/t)t>o on R''. The spatially homogeneous Landau (or Fokker-Planck-Landau) equation 
corresponds, in dimension d > 2, to the case where for some k : M_|_ i—t M_|_, 

(2) a^j{z) = KdzHdzpJy - z,Zj). 

Physically, one assumes that ^(r) = r'^^^, for some 7 G [—3, 1]. One talks of soft potentials when 
7 < 0, Maxwell molecules when 7 = 0, and hard potentials when 7 > 0. We consider in this paper 
the case of Maxwell molecules, or of pseudo-Maxwell molecules, where k is supposed to be smooth 
and bounded. 

This equation arises as a limit of the Boltzmann equation when all the collisions become grazing. 
Wc refer to Villani [TUIIIIIIII] and the many references therein for physical and mathematical details 
on this topic. See Cordier-Mancini and Buet-Cordier-Filbet Ij for a review on deterministic 
numerical methods to solve ([T]). 

1.2. Notation. Let V ~ 'P{M.'^) be the set of probability measures on Mf^, and Vk = {fJ- G 
V,mk{n) < 00}, where nikif^) = J \x\''n{dx). 

For x,y € M.'^, we set = i^fx^)^, and {x,y) = x*y = Y^fxtyi. We consider the norm 

= sup{|(Mx, x)!, \x\ = 1} = max{|A|, A eighenvalue of M} on Sd- Recall that for A e 
mf{{Ax, x), \x\ = 1} = l/|yl^^|. All A G admits a unique square root A'^ g S'^ , and we have 
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Definition 1. Consider a : R'^ i-^ S'j'. Let b : i-^ M.'^ be defined by bi{x) = Yl'j=i^j^ij{^)- 
Assume that \a{x)\ + |6(a;)| < C(l + (which is the case when a is defined by (0j with k e 

C^). A measurable family {Pt)t>o C V2 is said to be a weak solution to flj) if for all t > 0, 
sup[Q_t] 7712 (Pg) < 00 and for all ip £ C^(E'*), 

(3) / ^{x)Pt{dx)= I ip{x)Poidx)+ f ds f f Ps{dx)P,{dy)Lp{x,y), 

where Lip{x,y) = \ Y.tj=i O'iji^ " y)dfjf{x) + Y.t=i ~ y)diip{x). 

All the terms make sense due to our conditions on a, b, Pt. See Villani [11] for a similar formulation. 

1.3. Known results. To our knowledge, the first (and only) paper proving a rate of convergence 
for a numerical scheme to solve ([T]) is that of Fontbona-Guerin-Meleard [1]. Their method relies 
on a stochastic particle system. The aim of this paper is to go further in this direction. 
Let us thus recall briefly the method of |4Jj relying on the probabilistic interpretation of ([1]) 
developped by Funaki [5], Guerin [7]. 

Let a ■.R'^ 1-^ and 6 : R'^ i-^- R'' be Lipschitz continuous functions, and let Po G 7^2- A M^'-valued 
process (Xt)t>o is said to solve EQ{PQ,cr, b) if C{Xq) — Pq, and if for all i > 0, setting P^ = £{Xt), 



(4) Xt^Xo+ / a{Xs-x)Wp{dx,ds)+ / b{Xs ~ x)Ps{dx)ds. 

Jo Jr^ Jo Jmi^ 

Here Wp{dx,dt) is a M'^-valued white noise on [0, 00) x M.'^, independent of Xo, with independent 
coordinates, each of which having covariance measure Pt{dx)dt (see Walsh [14]). 
Existence and uniqueness in law for Eq(Pq, a,b) have been proved in Guerin [7]. If furthermore 
a{x)a*{x) — a{x) and bi{x) = J2j=i ^i'^y (^)j then {Pt)t>o is a weak solution to ([T]). The condition 
that (T and b are Lipschitz continuous is satisfied in the case of the Landau equation for Maxwell 
or pseudo-Maxwell molecules. 

In [3], one considers an exchangeable stochastic particle system {Xl'"')t>o,i=i,...,n, satisfying a 
S.D.E. driven by Brownian motions. It is then shown that one may find a coupling between a 
solution {Xf)t>o to Eo{Po, o", b) and such a particle system in such a way that 



E 



snp\Xl'--Xl\' 



< Crn 



-2/(d+4) 



under the condition that Pq has a finite moment of order d+b. The proof relies on a clever coupling 
between the the white noise and n Brownian motions. In particular, one has to assume that Pt 
has a density for all i > 0, in order to guarantee the uniqueness of some optimal couplings. 

1.4. Another approach. For a : R'^ S'j", & : R'* M'' satisfying |a(a;)| + \b{x)\ < C(l + |a;p) 
and n e 7'2(R'^), we introduce 

a{x,fi)= a{x-y)n{dy), b{x,ti)= b{x - y)fj.{dy). 

JW JRd 

For each x G M.'^, /i € V2, o,{x,ii) is a nonnegative symmetric matrix and thus admits an unique 
symmetric nonnegative square root 02 (a;, /i) := [a(a;,/i)]2. 

Denote by W^j the law of the d-dimensional Brownian motion, consider Pq € P2, and let (Xq, B) ~ 
Po ® Wrf. We say that a R'^-valued process {Xt)t>o solves Ei (Pq, a, b) (or Ei (Pq, a, 6, ATq, B) when 
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needed) if E[sup[o,T] l^tP] < °° for all T and if for all t > 0, setting Pt = C{Xt), 

(5) Xt^Xo+ f a-^{Xs,Ps)dBs+ f 5(X„P,)ds. 

Jo Jo 

This equation is nonlinear in the sense that its coefficients involve the law of the solution. Compared 
to (j4|), equation ([5]) is simpler, since it is driven by a finite-dimensional Brownian motion, and since 
the nonlinearity does not involve the driving process. However, one may check that at least formally, 
solutions to ([U and (O have the same law. The link with ([T]) relies on a simple application of the 
Ito formula. 

Remark 2. Let (Xt)t>o solve Ei{Po,a,b). Assume thathi and that \a{x)\ + \b{x)\ < 

C(l -|- \x\'^). Then {Pt)t>o i^{Xt))t>o is a weak solution to (Qp. 

The natural linearization of ([5]) consists of considering n particles (^i'")t>o.i=i, 

(6) xr - ^0 + «^ [xr. \ E ^x>i^ + b (^xr, i 

Here (Xq, are i.i.d. with law Pq (X) W^j. We thus use n Brownian motions. When 

linearizing ([5]), one needs to use Brownian motions, since the white noise is infinite dimensional. 
However, one may check that the solution to ^ and the particle system built in [3] have the same 
distribution (provided aa* = a in [4, Equation (4)]). 

1.5. Main results. The main result of this paper is the following. 

Theorem 3. Assume that b is Lipschitz continuous, that a is of class , with all its derivatives 
of order 2 bounded, and that Pq e 7^2 ■ 

(i) There is strong existence and uniqueness for Ei(Pq, a, b): for any {Xq, B) ^ Pq W^, there is 
an unique solution {Xt)t>o to Ei{PQ,a,b, Xq, B). 

(ii) Let (Xq, be i.i.d. with law PQ®'Wd. There is an unique solution {Xl'^)t>Q,i=i,...,n 
to IBj). Assume that Pq £ Vn, and consider the unique solution {Xl)t>Q to Ei{Po,a,b, Xq, B^). 
There is a constant Ct depending only on d, Pq, a,b,T such that 




(7) E 



snp\Xl^--XU' 

[0,T] 



^(^T min^n-i/^n-1 sup {1 + \a{x, Pt)-^])^ dt < Ctu-^/^ . 



In the general case, we thus prove a rate of convergence in which is faster than n^^^'^'^^^\ 

If we have some information on the nondegeneracy of a{x,Pt), then a^[x,^) is smooth around 
^, Ps, and we can get a better rate of convergence. 

Assume for example that a is uniformly elliptic (which is unfortunately not the case of since 
a{x)x = for all x £ R''). Then sup^ \cl{x, Pt)~'^\ < supj, |a(j/)^-'^| < oo, and we get a convergence 
rate in n~^. 

In the case of the Landau equation for true Maxwell molecules, we obtain the following result. 

Corollary 4. Consider the Landau equation for Maxwell molecules, where a is given by (0) with 
K = 1 and bi{x) = J2j=i ^j'^ij(^) ~ ^ Then a, b satisfy the assumptions of Theorem\^ 

Let Pq £Vi, and adopt the notation of Theorem\^ (ii) . 

(i) We have E[sup[o^j.] \X]'" - X^p] < CTW^il + \ogn). 

(ii) Set Xq — J xPQ{dx). If a{xQ, Pq) is invertible, then EfsupjQ \Xl'" — X^\'^] < Crn"^ . 
We finally consider the case of pseudo-Maxwell molecules. 
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Corollary 5. Consider the Landau equation for pseudo- Maxwell molecules, where a is given by 0j 
with K £ C'^(M-(-), and bi{x) = X]j=i ^^j^iji^) — ~('^ ~ l)K(|a;p)a;i . Assume that k' has a bounded 
support. Then a,b satisfy the assumptions of Theorem\^ (ii) . 

Assume furthermore that Pq S ^4 has a density with a finite entropy J PQ{x)logPo{x)dx < oo, 
and that k is bounded below by a positive constant. With the notation of Theorem we have 
E[sup[o,r]l^/'"-^t'l']<CTn-i. 

1.6. Time discretization. To get a simulable particle system, it remains to discretize time in 
Let > 1, and consider pn{s) = ^k>o ■^lse[fc/jv,(fc+i)/jv)- Consider the simulable particle 
system (Xt'"'^)t>o_i=i,...^„ defined by 



(8) xi: 



X' 



a3 X 



dBl 



b X^'^'f,, 



ds. 



Theorem 6. Assume that b is Lipschitz continuous, that a is of class , with all its derivatives 
of order 2 bounded, and that Pq £ T'2- Let {Xq, -B*)i=i....,n be i.i.d. with law Pq^'SM d- Consider the 
unique solutions {Xl'^)t>o,i=i,...,n to ^ and (^t'"'^)t>o,i=i,...,n to (0). Then there is a constant 
Ct depending only on d, Pq, a,b,T such that 



(9) 



E 



I -i^l.n vl.n.N \2 

[0,T] 



< CtN-^. 



1.7. Conclusion. Choosing for example a,b, and Pq eis in Corollary [J]- (ii) or as in Corollary [51 
denoting by {Pt)t>o = i^i^t))t>o the weak solution to the corresponding Landau equation, we 
obtain for any ip £ Cl,hy exchangeability. 



sup E 

[0,T] 



1 " r 

-Y^^iXr'")- / ^{x)P,idx 



Thus if one simulates the discretized particle system ([8]), and if one computes ^ J2i 'P(^t'"'^); 
get an approximation of J ip{x)Pt{dx), with a reasonnable error. 



1.8. Plan of the paper. In Section [2l we give the proofs of Theorems [3] and [G] Section [3] is 
devoted to the proofs of Corollaries 2] and El 

In Section 21 we briefly deal with the case of soft potentials, but our theoritical results do not 
extend well. Numerical results are given in Section [5l Finally an appendix lies at the end of the 
paper. 



2. General proofs 

In the whole section, we assume that Pq £ V2, that a : R'^ 1-^ is of class C^, with bounded 
derivatives of order two, and that & : M'' 1-^ R'^ is Lipschitz continuous. We denote by C (resp. 
Ct, Ct.p) a constant which depend only on a, 5, d, Pq (resp. additionally on T, on T,p) and whose 
value may change from line to line. 

For iJL,v £ V2, we set j/) = min {E [\X - Fp] ; C{X) = ^, C{Y) = 1^}. See Villani [T^ for 

many informations on the Wasserstein distance 'W2- 
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2.1. Preliminaries. Our results are mainly based on the two following Lemmas. 
Lemma 7. For all E V2, all x,y E M.'^, 

\J {x, ^l) - (y, + \b(x, /i) - b{v, < Ci\x - + WHii, v)), 

\a^{x,^l)\^ + \b{x,^i)\^ < C(l + m2(/i) + \xf). 

Proof. Step 1. For ji E V2 fixed, we consider the map A : R"* 5'j~ defined by A{x) — a{x,^). 
Then D'^A{x) — J^a D^a{x — y)ii{dy), is clearly uniformly bounded. Lemma flOl ensures us that 
||-D(A^)||oo is uniformly bounded, so that \a^{x,p) — a^{y,^)\ = \A^{x) — A^{y)\ < C\x — y\. 

Step 2. We now fix a; e W^, and consider iJi^v E 7^2- We introduce a couple {X,Y) of random 
variables such that C{X) = /i, C{Y) ^ v, and W|(/i, v) = E[|X - YW We define ^ : M 5*+ by 
Ait) = E [a{x - [tX + (1 - t)Y])]. Then A{{)) = ¥\a{x - Y)] = a{x, v) while ^(1) = E[a(.x - X)] = 
a{x,fj,). Furthermore, 

\D^A{t)\ = mix - YfD^aix - [tX + (1 - t)Y])]\ < \\D^a\\^M[\X - Y\^] ^ CWH^i, v). 
Lemma [TOl ensures us that ||(A^)'||oo < CyV2{f^,i^), so that \ai{x,fi) — ai{x,i^)\ — — 

A^o)\ < cmii^.v). 

Step 3. The growth estimate (for a) follows from the Lipschitz estimate, since |a^(0, (5o)P — 
|a2(0)p < 00, and since yV|(/x, (5o) = m2{pi)- 

Step 4- The case of b is much simpler. For /i, j/ £ V2, we introduce X,Y as in Step 2. Then 
\bix,^l)-biy,u)\^ = \E[b{x-X)-b{y-Y)]\^ < C{\x-y\^ +E[\X -Y\]^) < + t.)). 



The growth estimate follows from the Lipschitz estimate, since |6(0, <5o)p = |^(0)P < 



□ 



Lemma 8. Let Yi be i.i.d. 

E 



-valued random variables with common law fi E Vi- Then 



< C 



1 + mi{^) 



Proof. We denote by Ei the expectation concerning only Yi, and by E2.„ the expectation concerning 
only Y2, . . . , Yn. We observe that for alH = 2, . . . , n, we have a{Yi, fi) = E2^nla{Yi — Yi)], whence 
aiYi,fi) = E2,n[j^E2aiYi - Yi)]. We also have a{Yu ^r.) = ^Ei a(^i - Yi). As a 

consequence, 



E 



2Ei <^ E2, 



a(ri,^) 



< 2E 



^ n ^ n 

-Y^a{Y,-Y,) Y.a{Y,-Yi] 



^ n " 1 ^ 

— 1 71—1 



2-/^n ~f" 2t/j2 



An immediate computation, using that \a{x)\ < C(l + \x\^), shows that J„ < C(l + mi{y))/r 
On the other hand, since the random variables Yi — Yi are i.i.d. under E2,n, 

Jn <IEi < [^^Y.aki{Yi - YM \ < ^^Ei J J2 yar2,naki{Yi - Y2) 



<- 



C 



c 



c . 



-El {E2,„ [|a(ri - r2)|2] } < -E [|a(yi - ^2)1'] < -(1 + m4(M)), 
1 ^ '■ ^ ' n n 
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again since |a(a^)| < C(l + \x\'^)- The same computation holds for 6, replacing everywhere m4{^) 
by m2(Ai), since |6(a;)| < C(l + □ 

2.2. Convergence proofs. We start this subsection with some moment estimates. 

Lemma 9. (i) Let {Xt)t>o solve Ei{PQ,a,b). Assume that mp(Po) < for some p > 2. Then 
E[sup[o \Xt\P] < oo for all T > 0. 

(li) Let {Xl'")t>o,t=i,...,n solve 0j. For allO < s <t <T, E[\Xl'" ~ < Crlt - s\. 

Proof. Point (i). Set Pt=£{Xt). Using the Burkholder-Davies-Gundy inequality for the Brownian 
part, and the Holder inequality for the drift part, we obtain, for all < t < T, 



E 



snp\X,\P 

lo,t] 



<CpE[\Xo\P] + Cp dsE \a-2iX,,P,)\P + Cp^r dsE P,)n . 



But Lemma [7] implies that E[\a^Xs, Ps)]^ + |5(X,,P^)|p] < CpE[l + \Xs\p + 7712(^5)^/^]- Fur- 
thermore, since Ps = C{Xs) and p > 2, we deduce that m2{PsY^^ < E[|Xs|p]. As a conclusion, 
E[sup[o,(] \XsY'] < CpE[\Xq\p] + Cp^T /o dsE[l + \Xs\p], whence the result by the Gronwall Lemma. 

Point (a). Using the Cauchy-Scharz and Doob inequalities, we see that for < s < t < T, 



E 



<C / du. 



\aHxi^-,-Y,5^..)\' 



+ Ct du. 



(10) 



<Ct / duE 



(^E'^x-) <CTfduE[l + \Xl--\'] 



□ 



We used Lemma[7]and that E[m2(i X^i '^xi'")] = ^ Ei ]E[|A:;;"P] = IE[|A:1'"|2] by exchangeability. 
Applying (Uni) with s = 0, we get E[\Xl'"\^] < CE[\X^\'^]+Ct /q du[l+E[\Xl-"\^]du. The GronwaU 
Lemma allows us to conclude that supjQ ATj"^'"]^] < Ct- Applying a second time pO|). we 
deduce that EilXl'"" - X^-^l^] < Crlt - s\. 

Proof of Theorem\^ We consider Pg ^ P2 fixed. 
Point (i). Let (ATq, P) - Pq W<j. 

Uniqueness. Assume that we have two solutions X^Y to Pi(Po, a, 6, Xq, P), and set Pt — C{Xt) 
Qt — C{Yt). Using the Cauchy-Schwarz and Doob inequalities, we obtain, for < < < T, 



E 



SUp|X, 

[0,t] 



(11) 



<Ct / E[\a-2{Xs,Ps)~a^{YsMs)\'' + HXs,Ps)-b{Ys,Qs)\'']ds 
Jo 

<Ct I E [\Xs - + Wl{P,,Qs)\ ds<CT f E [\Xs ~ Y,f] ds. 
Jo Jo 



We used Lemma[7]and the obvious inequality yV|(P<;, Qg) < E[|Ars — Kjp]- The Gronwall Lemma 
allows us to conclude that X = Y. 

Existence. We consider the following Picard iteration: set X^ = Xq, and define, for rt > 0, t > 0, 



(12) 



xo+ aHx:,c{x:))dB,+ b{x:,c{x:))ds. 



We get as in (HU, for < i < T, E[sup[o,t] \X^+^ - AT^f ] < Ct j^EWX^ - Xl'-^\'^]ds. Thus 
there classically exists {Xt)t>o such that lim„E[sup[Q jAT" — Xt\'^] = for all T, which implies 
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that lim„ sup[Q yV|(£(X"), £(Xf )) = 0. Passing to the Umit in we see that X solves 

Ei{Po,a,b,Xo,B). 

Point (a). First of all, the strong existence and uniqueness for ^ follows from standard theory 
(see e.g. Stroock-Varadhan ^8 ), since for each i, the maps [xi, . . . ^Xn) i— > b{x^ , j^Y^ ^x'') and 
(xi , . . . , x„) (x% ^ Jj-fc ) are Lipschitz continuous (use Lemmas [7] and [T^]) . 

We now consider {Xq,B'^) i.i.d. with law Pq <8i W^, the solution (^f '")t>o,i=i,...,n to and 
for each i = l,...,n, the unique solution (X()t>o to Ei{PQ,a,b,XQ,B^). For each t > 0, let 
Pt — C{Xl) = ■ ■ ■ = C{X"). Due to the Cauchy-Schwarz and Doob inequalities, for < t < T, 



snp\Xl'--Xl\^ 
[o,t] 



dsl 



< Cj 



ds E 



A«(s) 



where 

A„(s) :=E 



(13) ^■.Al{s) + Alis). 

Using Lemmas [7l and fT2l we obtain, for < t < T, 



E 



sup l^i-" - X 

[0,t] 



1|2 



< Ct 



< 



ds 



(e ^xl--xl\^ + wl(^l^p,,.,'-±s,^ 



+ A„(s) 



^1 

< Ct / rfsE [IX]'" - Xlf] +Ct I dsAnis) 
Jo Jo 



Ct f rfsA„(s) 
Jo 



by exchangeability. The Gronwall Lemma ensures us that 



(14) 



E 



snp\Xl'^-Xl\' 

[0,T] 



<Ct [ dsAnis). 
Jo 



It remains to estimate A„(s). The random variables X^ , . . . , X" are i.i.d. with law Pg. Thus 
Lemma m shows that A^(s) < C(l + ■mi{Ps))/n < Ct/ti for s <T, due to Lemma|9]-(i) and since 
Pq G by assumption. Next, we use Lemma llll -fi). the Cauchy-Schwarz inequality, and then 
Lemma m for s < T, 



Alis) <E 



Xl-j2^xi]-a{Xl,C{Xl)) 



< C 



l + mi{Ps 



< 



Ct 
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But one may also use Lemma [TTI-fii) instead of Lemma [TTI-(i'). and this gives, for s <T, 



<Csup|a(:r,P,)-i| fii^^^) < ^sup|a(x,P.)-i| 



Thus A„(s) < Cxn ^ + Ct min(n ^^^,n ^ sup^ \a{x,Ps) ^|). Inserting this into (fT4)l . we obtain 
0. □ 

Proof of Theorem\^ Using Lemmas [71 and [T^ we get as usual (see (jlll) ). by exchangeability, 



E 



sup|Xl'"-Xi'"-^P 

[0,*] 



< Ct / E 



|^l,„_^l,nJV|2 
I s Pn[s) I 



< Ct / E 




3 



|^l,„_^l,nJV|2 
' " Piv(s)l „ 



1 " 



ds 



I s Pn(s) I 



<Ct E [l^l'" - + Ct 







|^l,«_^l,n |2 



ds. 

Using finally Lemma [9]-(ii), and since |s — pn{s)\ < 1/^, we deduce that E[|X]'" - ^p^",)^] < 
Ct/N. The Gronwall Lemma allows us to conclude. □ 

3. Ellipticity estimates 

We start with the 

Proof of Corollary^ Recall here that a is given by (O with k = 1 and b{z) = —{d — l)z. Thus h 
is Lipschitz continuous, and the second derivatives of a are clearly bounded. We consider a weak 
solution {Pt)t>Q to III). 

Simple computations using ([3]) (with (p{x) = Xi, ip{x) = jxp) show that dt J xPt{dx) = and 
dtm2{Pt) = 0. We classically may assume without loss of generality that / xPt{dx) — J xPo{dx) — 
0. We also assume that ■m2{Pt) = m2{Po) > (else — X^'" — a.s.). 

We now bound from below {a{x, Pt)y, y) for x, y e R'', t>Q. 

A simple computation, using that / xPt{dx) = 0, shows that a{x,Pt) — a{x) + a(0,P(). Thus for 
all t > 0, x,y G K'', setting m2 (Pt) = J XiXjPt{dx) 

{a{x,Pt)y,y) > {a{0,Pt)y,y) = ^ y.y^- [7712(^4)% - (Pt)] = "^2(^o)|yP -^yty^rn^ {Pt). 



Using ([3]) with Lp{x) = XiXj, we deduce that 

dtmiiPt) = 2m2{Pt)5^, - 2dfT^^ {Pt) = 2m2(Po)% - 2dm^^(Pt). 

We thus obtain 

dt{a{Q,Pt)y,y) = -Y^y^y^dtm^ {Pt) = 2dY,y,y,m!-i {Pt) - 2m2(Po)|y|' 

i,j hi 

^2{d - l)m2(Po)|2/|2 - 2d(a(0, P*)?/, 2/). 
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Set Ao = inf{(a(0,Po)y,2/), \y\ = 1} > and Ai = ^m2(Po) > 0. For all t > 0, all x,y e M'', 

{a{x, Pt)y, y) >(a(0, Pt)y, y) = (a(0, Po)y, y)e-^^' + \i\yWl - e'^"') 
(15) >|ynAoe-2''* + Ai(l-e-2'^*)]. 

We now prove point (i). We deduce from that 

(a(a;,P,)2/,y)> Ai(l-e-2'i*)|y|2. 

As a consequence, \a{x,Pt)~'^\ < l/[Ai(l — e^^''*)] < C/t + C. Inserting this into ([7]), we get 
E[sup[o,T] l^t^ -^^"P] < Ct /o^niin(7i-i/2^„-i + {nt)-^)dt < CTn-i(l + logn). 

To get (ii), we use (fTSl) and that by assumption, Aq > 0. We deduce that 

{a{x,Pt)y,y) > lyl^Xoe^^''' + Ai(l - e^^d*)) > ^^^^^^^ \^)\y\y2. 
As a consequence, \a{x, Pt)^^\ < 2/ min(Ao, Ai). Inserting this into we get E[sup[o.T] ~ 
Xl'y] < Ct Jo mm{n-^/\n-^)dt < Ctu-K □ 
It remains to give the 

Proof of Corollary \^ Recall here that aij(x) — K(|a;p)(|xp(5y — XiXj) and that h{x) = —{d — 
that K is and that k' has a bounded support, so that a has bounded derivatives of 
order 2, and h is Lipschitz continuous. We consider a weak solution {Pt)t>o to ((J). As previously, 
we classically have m2{Pt) = m2(Po)- Furthermore, it is again classical and widely used that the 
entropy of P* is non-increasing, so that J Pt{x) \ogPt{x)dx < J Po{x)logPo{x)dx = C < oo for all 
times, see Villani [lOl HH HI] . 

If we prove that there is Aq > such that for all t > 0, x,y & R'', (a(x, Pt)y, y) > Ao|yp, then we 
deduce that \a{x,Pt)~^\ is uniformly bounded, so that the Corollary follows from 
Observe that setting aij{x) = \x\'^5ij — XiXj, we have {a{x, Pt)y,y) > Xi{a{x, Pt)y,y), where 
Ai > is a lowerbound of k. But it is shown in Desvillettes- Villani [31 Proposition 4] that for 
Eq G M.^,Hq G K_|_, there is a constant ceo,Ho > such that for any probability density function 
/ on R''- such that m2(/) < Eq and / f{x)\og f{x)dx < Hq, {a{x,f)y,y) > CEo,Ho\y\'^- Actually, 
they consider the case where aij{x) — \xp {\x\'^Sij — XiXj) for some 7 > 0, but one can check that 
their proof works without modification when 7 = 0. We finally obtain (a(x, Pt)y, y) > XiCEo.Ho\y\'^ 
for all t > 0, x,y G M.'^ , which concludes the proof. □ 

4. On soft POTENTIALS 

We consider in this section the spatially homogeneous Landau equation for soft potentials, which 
writes ([1]) with aij{z) = \z\''{\z\'^6ij—ZiZj) for some 7 € [—3, 0), the Coulomb case 7 — —3 being the 
most interesting from a physical point of view. Then we have bi{z) = djOij^z) = —{d— l)\zp Zi. 

Simulation with cutoff. We restrict our study to the case where 7 G (—2, 0]. We assume that Pq 
has finite moments of all orders, and has a density with a finite entropy / Po{x) \ogPo{x)dx < 00. 
For £ > let Ke : M+ 1-^ M+ of class C^, nondecreasing, with Ke{z) = z for z > e, Ke{z) = e/2 for 
z g [0,e/2], with +e|K"(z)| < C. Consider then a^,bs be defined as a, 6 with \z\'' replaced 

by [Kgdzj)]'''. Then is of class C^, with all its derivatives of order 2 bounded by Ce^, and is 
Lipschitz continuous with Lipschitz constant Ce'^ . 

We thus may apply Corollary O and Theorem [S) Denote by {P^)t>o ~ {C{Xl'^))t>o a weak 
solution to ^ with and Pq = Pq. Then we believe that our results, plus some moment and 
ellipticity estimates (uniform in e e (0, 1]), will give something like E[sup[Q ^x^'^'^'^^ — Xl''^\'^] < 
(n-i + Af-i)exp(CTe2^), where (X*'"^"^^") t>0A=i,..,n solves ([5]) with a^,b^ instead of a, 6. 
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On the other hand, we may apply the techniques introduced in [S] to estimate yV|(Pt, ), where 
{Pt)t>o = {^{^t))t>o is a weak solution to ([T]) with a and Pq. We believe that, with a convenient 
coupling, it is possible to obtain something like supjQ j.^ ¥,[1x1'^ — X^l'^] < Ct£^- 

One would thus get sup^o Ellxl'""'^'' - X}\^] < C't (e^ + {n-^ + N-'^) e^^^^^Y This is of course 
an awfuU rate of convergence. It does not seem reasonnable to handle a rigorous proof. 

Simulation vifithout cutoff. However, the particle system ([S]) is still well-defined and simulable 
for soft potentials (with 7 G [—3, 0]), at least if we replace ^ S^^k.n.N by ^ X^fe^^i S^k.n.N and if Pq 
has a density. Based on the well-posedness result of [5], we hope that, at least when 76 (—2, 0], one 
might obtain the same estimates as in Corollary [5] and Theorem IH] ( under additionnal conditions 
on Pq). The proof however seems to be quite difhcult: we do not know how to get a sufficiently 
good estimate of quantities like jX^'"'^ — 

5. Numerics 

Let us first observe that for the Landau equation [1] where a is given by ^ and bi = J^j ^jCiij, the 
simulable particle system ([5]) is conservative, in the sense that it preserves, in mean, momentum and 
kinetic energy: for alH = 1, . . . , n, alH > 0, E[X('"'^] — J xPo{dx) and E[|Xf '"'^p] — m2(Po). 

We consider here the Landau equation for soft potentials, for some 7 G [—3,0], described in the 
previous section, in dimension d = 2. We use no cutoff procedure in the case 7 < 0. We consider 
the initial condition Pq with density Po{xi, X2) = f{xi)g{x2), where / is the Gaussian density with 
mean and variance 0.1, while g{x) = {f{x — 1) + f{x + l))/2. The momentum and energy of Pq 
are given by (0,0) and 1.02. 

Thus in large time, the solution Pj should converge to the Gaussian distribution with mean (0, 0) 
and covariance matrix O.5I/2, see Villani [T^]. 

We use the particle system (jS]) with n particles, and N steps per unit of time. Easy considerations 
show that the computation of ([8]) until time T is essentially proportionnal to TNn? , and should not 
depend too much on 7. However, it is consequently faster when 7 = for obvious computational 
reasons. Let us also remark that the law of ([S]) does not change when replacing a 2 by any a such 
that a(x, ij)a{x, jj)* = a(x,ij). We thus use a Cholesky decomposition, which is numerically quite 
fast. Let us give an idea of the time needed to perform one time-step: with 7 = 0, it takes around 
7.10-3 seconds (n = 500), 0.15 s (n = 2500), 3.5 s {n = 12500), and 13 s (n = 25000). The 
computations are around 10 times slower when 7 < 0. 

Now we alway use n = 5000 particles, and N — 200 steps per unit of time. We draw, for different 
values oft and 7, the histogram (with 80 sticks) based on the second coordinates of (X('"'^)i=i,...^„. 
The plain curve is the expected asymptotic Gaussian density, with mean and variance 0.51. The 
convergence to equilibrium seems to be slower and slower as 7 is more and more negative. 



gamma=0, t=0.25 



gamma=0, t=0.5 



gamma=0, t=1 
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For too small values of 7 (say 7 < —2.5), the numerical results are not so convincing. This is not 
surprising, since the coefficients are more and more singular as 7 becomes smaller and smaller. 

6. Appendix 

The following Lemma can be found in Stroock-Varadhan (when p ~ d) 8, Theorem 5.2.3], or in 
Villani 9, Theorem 1] (for a more refined statement including all possible values of p and d). 

Lemma 10. Let A : M.^ i—> 5^, for some p > 1, d > 1, be of class , with all its derivatives of 
order 2 bounded. Then \ \D{A^)\\oo < Cp^d\/\ \D'^A\\ao, where Cp^d depends only on p and d. 

We also need the following estimates, which are probably standard. 

Lemma 11. For A, B e Sj, 

(i) there holds \ A^ -B^\< ^\A- B\ 

(li) and 1^3 ~B^< 0Tiin(|^-i|,|B-i|) x |^ - B\. 

Proof. We start with point (i). Let a = \A^ ~ B^\. There is a unit vector e € R'' such that 
1(^2 — i?2)e| = ere, and we may assume that (^2 — B2)e = ae (else, change the roles of A,B). 
Then, using that B^ is nonnegative, 

\A-B\> {Ae - Be, e) = {{A^ - B^)e, {A^ + B^)e) = (ere, ae + 2B^e) > a^\e\^ = . 

We now prove (ii). First observe that {A^x,x) > |xp/|v4~^/^| for all x E R'^. As previously, 

\A-B\ >iiAi ~ B^)e, {A^ + B^)e) = cr(A^e,e) + cr(B^e,e) 

>a|e| Vl^-^/'l + ^\e\V\B-^'^\ = a j ^f\A^\ + a j ^f\B=^\ > a/Vmin(|A-i|, 

whence \Ai - B^ ^ a < y^min{\A-^,\B-^)\A ~ B\. □ 

We conclude this annex with an elementary fact on the Wasserstein distance. 

Lemma 12. For xi, . . . , x„, yi, . . . , y„ e WHiEi^.^^Ei^v.) < 71 ^1 - y^? ■ 

Proof. Let U be uniformly distributed on {1, . . . , n}, set X — xjj a.n Y = yu . Then A ~ — , 
El K > and E[|A - Y\^\ = i Ei " " □ 
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